Scaling Ion Trap Quantum Computation through Fast Quantum Gates 
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We propose a method to achieve scalable quantum com- 
putation based on fast quantum gates on an array of trapped 
ions, without the requirement of ion shuttling. Conditional 
quantum gates are obtained for any neighboring ions through 
spin-dependent acceleration of the ions from periodic photon 
kicks. The gates are shown to be robust to influence of all the 
other ions in the array and insensitive to the ions' tempera- 
ture. 

PACS numbers: 03.67.Lx, 03.65.Vf, 03.67.Pp, 32.60.Qk 

Trapped ions constitute one of the most promising sys- 
tems for implementation of quantum computation. Sig- 
nificant theoretical and experimental advances have been 
achieved for this system [1-8], and quantum gates have 
been demonstrated at the level of a few ions [5-8]. The 
current central problem is to find methods to scale up 
this system for larger-scale quantum computation [9-14] . 
A particularly promising approach to scalability is based 
on shuttling ions in complex traps to different regions for 
storage and for quantum gate operations [9-11]. Interest- 
ing initial experiments have been reported on separating, 
shuttling, and sympathetic cooling of the ions [15-17]. 
Nevertheless, these experiments also indicate that fast 
separation of the target ions is a challenging task [15], 
which limits the speed of any collective quantum gate in 
a scalable structure. 

In this paper, we propose a scaling method based on 
fast quantum gates on an array of trapped ions. Very 
recently, Garcia-Ripoll, Zoller, and Cirac proposed a re- 
markable scheme of two-bit quantum gates [18], which 
can operator faster than the ion trap frequency, where the 
latter was thought before as a speed limit of the ion gates. 
The original fast gate in [18] was designed for two ions, 
with in mind that one achieves scalability through the ion 
shuttling. Here, we propose an efficient scaling method 
for ion trap quantum computation based on the concept 
of fast quantum gates. This scaling method provides an 
alternative scheme to achieve scalability, avoiding a series 
of challenges associated with the ion shuttling. 

The extension of fast quantum gates from two ions to 
a large array of ions is actually quite challenging: even if 
the laser pulses are shined only on two ions, all the ions 
in the array affect each other through the long-range Col- 
umn interactions, and all the phonon modes need to be 
taken into account as the motional sideband addressing 
is not possible with a fast gate. The mutual strong influ- 
ence between the ions normally sets a significant obsta- 
cle to scalable quantum computation. However, we show 
that as long as the gate speed is faster than the local 
ion oscillation frequency (specified below) , this unwanted 
influence can be arbitrarily reduced with a remarkable 
method for noise cancellation. Besides the proposal of an 



efficient scaling method, we also give a different design 
of the fast quantum gates with the use of only periodic 
laser pulses. Such periodic pulses are significantly easier 
for experimental realization. 

We consider an array of trapped ions, which could be in 
any convenient geometry. For instance, Fig. la shows a 
possible configuration of the array from multi-connected 
linear Paul traps. The qubit is represented by two sta- 
tionary internal states of the ion, denoted as |0) and |1) 
in general. We assume that the distance d between the 
neighboring ions is appreciably large ( d > 5/Ltm) so that 
single-bit operations and measurements can be done with 
no difficulty. The critical issue is then to implement pro- 
totypical two-bit gates for any neighboring ions, which, 
combined with simple single-bit gates, realize universal 
quantum computation. 
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FIG. 1. la. A possible geometry of the ion array for scal- 
able quantum computing from multi-connected linear Paul 
traps. Conditional quantum gates are achieved for any neigh- 
boring ions though application of a spin-dependent kicking 
force on these two ions. lb. The typical trajectory of the 
kicked ion in phase space, where "+" denotes the initial po- 
sition, "o" and "x" denote the final position of the ion in 
|0) or |1) states, which in general is different from the initial 
position as the ion has an initial velocity. 

For the purpose of the two-bit gate on any neighbor- 
ing ions, say, i and j, we assume that one can mechan- 
ically accelerate each ion, with the acceleration direc- 
tion depending on its internal state |0) or |1) (as illus- 
trated in Fig. la). So the acceleration force is writ- 
ten as F tJ = F^ a =i,j wh ere a' a = \l) a <1| - |0) Q (0| 
is the Pauli operator. The state-dependent acceleration 
can be realized, for instance, through coherent photon 
kicks from fast laser pulses [18]. A pair of short Raman 
pulses with the wave vectors kiand k2 applied on the 
ions i and j will give a state-dependent momentum kick 
h (ki — k 2 ) a z a (a = to the ions for each 7r-rotation 
of the states |0) and |1) . For convenience, we assume 
ki _L k 2 and |ki| s» | | = k c = 2tt/X c , where A c is the 
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carrier wave length. The kiand k 2 have a 45° angle to the 
two-ion axis and we alternatively apply the Raman pulses 
in the (ki,k 2 ) and (— ki, — k 2 ) directions. Each pulse se- 
quence is periodic with a repetition frequency f r /2 much 
larger than the ion trap frequency, so the net acceleration 
force on the ions is along the two-ion axis with the mean 
magnitude F = \2hk c f r . The momentum kicks from 
the periodic laser pulses are actually discrete in time, but 
their effect is well approximated by a continuous kicking 
force F if more than ten pulses are applied during each 
gate, as we have checked in our calculations. 

To characterize this quantum gate, we need to take 
into account the Column interactions between all the 
ions. The conventional approach is based on decompo- 
sition into the canonical phonon modes. However, for a 
large number of ions, this approach becomes very com- 
plicated even if one just wants to find out all the phonon 
modes. Here, we develop a different theoretical frame- 
work which is more appropriate for description of the 
fast gates on a large ion array. We note that with fast 
laser pulses, the acceleration force F^ can be significantly 
larger than the forces from the Column interactions and 
the external trapping potentials. We therefore write the 
total Hamiltonian H into two parts H = Hq+Hi, where 
the dominant part 



H = -F(t) J2 



a=t,j 
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(1) 



accounts for the kicking force F (t) on the target ions i,j 
and the kinetic energy summarized over all the ions k 
(including i,j); and the second part 



Hi = Vfefe' (|a; fc - Xfc/|) + ^ V T (x k ) 



(2) 



k,k> 



accounts for the Column interactions Vkk' (\xk ~ x k'\) be- 
tween every pair of ions k, k' and the external trapping 
potentials Vr (xk) on each ion k. In Eqs. (1) and (2), 
Xk denotes the displacement operator of the fcth ion from 
its equilibrium position (so d Xk Hi (xk) = by definition) 
and pk is the corresponding momentum operator [19]. 

To first entangle and then disentangle the ions' inter- 
nal and motional states, we reverse the direction of the 
kicking force by two times as depicted in Fig. 2a. So 
the force F (t) in H is time dependent, with F (t) = F 
for < t < ti and 3i x < t < 4*i, and F (t) = -F for 
ti < t < 3ii (see Fig. 2a, the value of t\ will be specified 
later). Under these laser kicks, Fig. lb shows the typi- 
cal state-dependent trajectory of the i or j ion in phase 
space under the "free" Hamiltonian Hq. Interestingly, 
independent of the initial position or velocity of the ion, 
the wave packets corresponding to the internal states |0) 
and |1) first split and then rejoin after a round trip. Un- 
der the Hamiltonian Hq, the position operator Xk of the 
fcth ion evolves as 



Xk (t) = x k + (t/m)pk 



*(*) 



E 



(3) 



where Xk and pk are the initial position and momentum 
operators, and s (t) = L J Q Tl [F (t 2 ) /m] c?T 2 dTi. 
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FIG. 2. The sequence of the spin dependent kicking force 
for the one-cycle (2a), two-cycle (2b) and four-cycle (2c) 
schemes, with t 2 = ti/2 1/5 and t 3 = ti/2 2/5 . The solid 
(dashed) curves are respectively for the ions in |0) or |1) states. 

To investigate the internal state evolution of the 
ions, we transfer to the interaction picture with re- 
spect to the "free" Hamiltonian Hq. The Hamilto- 
nian Hi in the interaction picture becomes Hi (t) — 
Ul (t) Hi ({x k }) U (t) = Hi ({x k («)}), where U (t) = 
r|exp —i J * H a (t) /hdr | (T {■ ■ ■} represents the 

time-ordered integration) and Xk (t) = Uq (t) XkUo (t) 
is expressed as Eq. (3). The evolution opera- 
tor in the interaction picture is given by Ui (t) — 



T jexp -i f Hi{ T )/fidT \. 



Note that at the end of 

the gate (t = 4ti), the Uq (t) becomes an identity opera- 
tors for the internal dynamics (as s (4£i) = in Eq. (3)), 
therefore the gate operation is wholly determined by the 
evolution operator E/j (t) at time t = Ati. 

To have an expression of C7j (t) , we note that there are 
three different length scales in our problem. They are the 
distance between the neighboring ions denoted as d, the 
length scale of the conditional displacement s (t) denoted 
as s, and the magnitudes of Xk and (t/m)pk which are 
estimated by the ion oscillation length and denoted as £. 
For fast quantum gates, typically we have d s £. 
Therefore, to the lowest (0th) order of the parameter 
£/s, the coordinates x a (t) (a = are approximated 
by s (t) <t„. In this limit, at time t = 4ti, the internal 
and external dynamics of Ui (4ti) become disentangled, 
and the internal state evolution is described by 



U m = exp 



41 1 



(t) dial a) 



(4) 



where the derivative of the Coulomb potential = 
d Xi d Xj Vij (\xi — Xj\) ~ 2hc/ (137c? 3 ) (c is the light veloc- 
ity). In writing Eq. (4), we have expanded the interac- 
tion Hamiltonian Hi (t) to the second order of the small 
parameter s/d (the harmonic approximation). Equation 
(4) represents an ideal controlled phase flip gate on the 
neighboring ions i,j if (V%/h) ,/ 1 s 2 (t) dt = tt/4, which 
determines the gate time 
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T g = Ah ~ 7.09 [f?v 2 r c/d 3 ] 



-1/5 



(5) 



In deriving Eq. (5), we have substituted the expression 
of the kicking force F into s (t), and v r — hk c /m denotes 
the atom recoil velocity. Figure (3a) shows the gate time 
T g as a function of the repetition frequency f r of the 
laser pulses for the 111 Cd + , 40 Ca + and 9 Be + ions. The 
gate speed (1/T 9 ) increases with the kicking frequency 

as fr^° and decreases with the neighboring-ion distance 
as <i -3 / 5 . The total number of the laser pulses for each 
gate is given by N = f r T g , which increases with the gate 
speed as (1/T 9 ) 3/2 . We can also expand the Hamiltonian 
Hj (t) to higher orders of the small parameter s/d (be- 
yond the harmonic approximation) , which give some tiny 
correction to the internal dynamics Ui„ [20]. However, as 
long as we are in the lowest order of the parameter £/s, 
the internal and external dynamics become disentangled 
after the pulse sequence, and there is no intrinsic noise 
to the gate operation. 
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FIG. 3. 3a. The gate time T g for the controlled phase flip 
operation is shown as a function of the repetition frequency 
of the kicking laser pulses for the 111 Cd + (solid curve), 40 Ca + 
(dotted curve), and 9 Be + (dashed curve) ions, where the 
neighboring ions' distance d is taken as d = 10 fim, and the 
carrier wave length for the nl Cd+, 40 Ca+, and 9 Be + ions are 
assumed to be 215 nm, 393 nm, and 313 nm, respectively. 2b. 
The scaling parameter uJiT g in the gate fidelity is shown as a 
function of the neighboring ions' distance d for 111 Cd + with 
the pulse repetition frequency f r = l GHz (dashed curve) and 
f r = 10 GHz (solid curve), respectively. The curve is well fit 
by the scaling uiiT g cx d' — 9/10). 

We now examine the influence of the ions' oscillations 
by expanding the interaction Hamiltonian Hi (t) to the 
next (1st) order of £/s, which gives some intrinsic noise 
to the above conditional quantum gate after the ions' 
motional states are traced. This intrinsic noise comes 
from two respects: (i) The internal and motional states 
of the i and j ions do not become fully disentangled after 
the pulse sequence; (ii) Due to the mutual Coulomb in- 
teractions, the state-dependent trajectory of the i,j ions 
yields a state-dependent force on all the other ions, which 
entangles the gate operation with the motional states of 
all the ions in the array. When the Hamiltonian Hj (t) 
is expanded to the next-order of £/s, the evolution oper- 
ator Ui (4ti) after the pulse sequence is decomposed as 
a product of the dominant part given by U{ n in Eq. (4), 
and the noise part U ns given by 



exp 



^ cr« {fiakXk + lakPk) 



(6) 



7a/c = 

tor H» ak 



The coefficients (5 ak = (1 - 8 ak /2) H'{ ak J^ 1 s (t) dt, 

(1 - S ak /2) H'{ ak ff 1 (t/m) s (i) dt, and the fac- 
d Xa d Xfi Hi denotes the potential derivative 
at the ions' equilibrium position. The operator U ns en- 
tangles the internal state of the target ions with the 
motional states of all the ions in the array, and intro- 
duces a gate infidelity after we take trace over the mo- 
tional states. To give a quantitative estimate of this in- 
trinsic gate infidelity, we take the following simplifica- 
tions: (i) We denote ujl = y/H" kk /m as the local os- 
cillation frequency of the ion k (with all the other ions 
fixed on their equilibrium positions) and assume that 
this frequency is roughly the same for all the ions {u>l 
is fc-independent); (ii) The states of the local oscillation 
modes (x k ,Pk) are assumed to be thermal with the mean 
phonon number estimated by n = ksTt/ (fiu>i), where T t 
is the temperature; (iii) The gate fidelity F g is defined as 
F g = (*o| U} n p r U in |^o), where |^ ) is the initial state 
of thei, j ions, taken to be (|0) i + (|0) 3 . + |1>^ /2 

as an example, and the final reduced internal state 



P 



Tr„ 



U na U in \* ) (*o\U} n Ul 



is obtained by tak- 
ing trace over the motional states; (iv) We define the 
parameters £,s quantitatively as £ 

4tl O r,\ „1 ; f f4ti 



Csm)dt i j^s(t)dt 



y/h/ (tbwl) and 
With these reason- 



able simplifications, we find that the gate fidelity F g 



is estimated by F g ps exp 
the dimensionless constant 



-ci(n + l/2)(g/s) ,where 
ci is defined as c\ = 



{H'>j2H>> i3 y + E Ml (H[' ik /H^Y\ /4. With the 
(5) for T g , we can also ex- 



definitions of £, s and Eq. 
press the gate fidelity as 



F g s» exp [-C2 (n + 1/2) ujlT 9 



(7) 



The dimensionless constant C2 — 0.83ci (H"ij / H'u-t) ■ To 
estimate the typical values of c\ and C2, in calculating 
H"^ we take into account only the contributions from the 
intrinsic Coulomb interactions (neglecting the contribu- 
tion from the external trapping potentials) , and find that 
ci w 8.6 (12.0) and c 2 ~ 3.0 (3.0) respectively for a 1 (2)- 
dimensional infinite ion array. Note that the values of C\ 
and C2 are well bounded from above even if we summarize 
over an infinite ion array, resulting form the fact that the 
noise comes from an effective mutual dipole interaction 
between the ions (described by the second derivatives of 
Hi), which falls off very rapidly with the ions' distance. 

The critical parameter io^Tg is shown in Fig. 3b as a 
function of the distance d between the neighboring ions 
with the laser repetition frequency f r = 1 (10) GHz. Al- 
though this parameter scales down with increase of the 
ions' distance, its value is not tiny in the typical param- 
eter region, and it is ineffective to reduce this parameter 
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by further increase of the laser repetition frequency due 

— 2/5 

to the slow scaling u>lT 9 ~ f r .A significant intrinsic 
infidelity for the quantum gate seems to be unavoidable 
if we substitute the typical value of u){F g into Eq. (7). 
Fortunately, this is not the case as there is an elegant 
way to greatly reduce this noise. In the above analy- 
ses, we applied one cycle of the kicking force (shown in 
Fig. 2a) which pushes the ions to the left (right) side 
if they are in the |0) (|1)) state (see Fig. 1). We can 
improve the scheme by using a two-cycle force as shown 
in Fig. 2b. Each cycle of the evolution contributes the 
same amount of conditional phase which accumulated for 
a conditional phase flip gate (so in Fig. 2 we should 
choose t 2 = ii/2 1 / 5 = Tg/4/2 - 2 ). However, for these 
two cycles, the ions are pushed to the reverse directions, 
and the coefficients fl a k and "f a k in the noise operator 
U ns (6) thus have a reverse sign. Due to this sign flip, 
the noise effects from these two cycles exactly cancel with 
each other! This noise cancellation is perfect up to the 
1st order of the parameter £/s. To estimate the residue 
noise with a two-cycle force, we take into account all 
the higher order contributions which basically rotate the 
operators Xk and pk in Eq. (6) with a rotational angle 
6 s=s 4cjl£ 2 ~ 0.87ujLTg for two different cycles, so after a 
partial cancellation of the noise, the accumulated coeffi- 
cients flak and 7„fe in the noise operator U ns are reduced 
effectively by a factor of 1 — e l6 w iO, and the final gate 
fidelity F g 2 for the two-cycle scheme thus has the form 

F g2 w cxp -0.87 2 c 2 (n + 1/2) (uj L T g f . We can further 

reduce the noise by using more cycles of the kicking force 
to get noise cancellation up to a higher order of the ro- 
tation angle 9. For instance, Fig. 3c shows a four-cycle 
kicking force, which reduces the accumulated noise coef- 
ficients (3 a k and -fak by a factor of 2 [u)lT 3 ) x 4 -2 / 5 , and 
the corresponding gate infidelity F g 4 is given by 



F, 



g4 « exp 



-1.32C2 (n+1/2) {u L T g y 



(8) 



With such a improvement, the gate infidelity SF = 
1 — F g 4 has been tiny in the typical parameter region. 
For instance, SF w 0.006% with the laser repetition fre- 
quency f r w 1 GHz, the neighboring ions' distance d w 50 
fim, and thermal phonon number n <~ 1. 

Summing up, we have proposed a scheme to achieve 
scalable quantum computation with trapped ions based 
on control from fast laser pulses. This scheme has the 
following distinctive features: (i) As the ion shuttling is 
not required in this scheme, one can use an ion array 
in any convenient geometry from any types of ion traps; 
(ii) The scheme is insensitive to the temperature of the 
ions, and requires no challenging cooling of the ion crys- 
tal to attain the Lamb-Dicke limit. The temperature T t 
affects the gate fidelity through the mean phonon number 
n = ksTt/ (?iuil) in Eq. (8), and this influence is pretty 
weak. For instance, even with a hundred of phonon ex- 
citations (n ~ 100) right after the Doppler cooling, the 
gate infidelity is still well below 1%. (iii) The conditional 



gates are very fast even if we have a pretty large distance 
between the neighboring ions which allows easy separate 
addressing; (iv) The whole computation is also fast as 
the slow step of separating the ions is not required for 
achieving scalability in this scheme. The basic experi- 
mental challenge is to well control the fast laser pulses 
to induce series of coherent photon kicks. Periodic laser 
pulses with a repetition frequency varying from hundreds 
of MHz to few THz have been reported in many exper- 
iments [21]. These impressive achievements and rapid 
progress in control of short laser pulses, combined with 
the ion trap technology, indicate realistic prospects for 
building scalable ion trap quantum computers based on 
this approach. 
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